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ABSTRACT 

Warm dark matter (WDM) and self-interacting dark matter (SIDM) are often motivated by the 
inferred cores in the dark matter halos of low surface brightness (LSB) galaxies. We test thermal 
WDM, non-thermal WDM, and SIDM using high-resolution rotation curves of nine LSB galaxies. 
We fit these dark matter models to the data and determine the halo core radii and central densities. 
While the minimum core size in WDM models is predicted to decrease with halo mass, we find that 
the inferred core radii increase with halo mass and also cannot be explained with a single value of the 
primordial phase space density. Moreover, if the core size is set by WDM particle properties, then 
even the smallest cores we infer would require primordial phase space density values that are orders 
of magnitude smaller than lower limits obtained from the Lyman alpha forest power spectra. We also 
find that the dark matter halo core densities vary by a factor of about 30 from system to system while 
showing no systematic trend with the maximum rotation velocity of the galaxy. This strongly argues 
against the core size being directly set by large self-interactions (scattering or annihilation) of dark 
matter. We therefore conclude that the inferred cores do not provide motivation to prefer WDM or 
SIDM over other dark matter models. 

Subject headings: cosmology: observations — cosmology: theory — dark matter — galaxies: kine- 
matics and dynamics 



1. INTRODUCTION 

In the prevailing theory of galaxy formation, galaxies 
assemble inside cuspy cold dark matter (CDM) halos. 
Dissipationless CDM simulations show that these 
halos have steeply rising central densities that roll 
over from p ~ r _1 at ~ 1 kpc scales to p ~ r -°- 8 a t 
~ 100 pc scales (e.g., INavarro et al J 120041 : iGraham et al.l 
120061 INavarro et all 120081 ). There has been contin- 
ued debate in the literature over this prediction as 
numerous observational results indicate that rotation 
curve data are often more consistent with dark matter 
halos having a roughl y constant density core (e.g., 
Flores fe Prhnackl H99l iMoord H99l iMcGaugh et all 



20011: iMarchesini et all 120021: iGentile et all 
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Simon et all 1 20051: iKuzio de Narav et all 120061 
de Blok et alll2008tlde"Bloldl2009l ). 

The mismatch between the dissipationless CDM sim- 
ulations and galaxy rotation curve data has motivated 
a serious exploration of warm da r k matter (WDM) 
mode ls (iHogan fe Dalcantonl I200C; lAvila-Reese et all 
20011 lAbazaiian et all 120011: iKaplinghatl 120051 
Cembranos et al.l 120051 : iStrigari et all 120071) and self- 
interacting dark matter (SID M) (jSpergel fc Steinhardtl 
120001 Kaplinghat et al.l 120001 ) as alternatives to CDM. 
The goal is to maintain the success of CDM on 
large scales while producing cores in the dark matter 
distribution of small halos. 

Two generic classes of WDM particles are thermal par- 
ticles that are light and which kinetically decouple when 
relativistic (|Blumenthal et al1ll982tlDodelson fc Widrowl 
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Il994t lAsaka et al.ll2006h . and particles that are as mas- 
sive as typical CDM particles but populated by decays 
in the early universe ((Kaplinghat 2005; C embranos et all 
2005). In both classes there are regions of parame- 
ter space that are endowed with large particle velocities 
and correspondingly low primordial phase space densi- 
ties Q p = p/cr 3 . The initial phase space density im- 
poses a limit on the central phase space density of col- 
lapsed WDM halos for both thermal WDM particles 
(|Tremaine fc Gunnlll979f) and non-therm al particles aris- 
ing from earlv-decavs (|Kaplinghatll2005D . Warmer dark 
matter models have lower Q p and correspondingly larger 
limiting core sizes at fixed halo mass. If the cores in 
dark matter halos are set by primordial phase-space con- 
straints, then for a given Q p there are predicted rela- 
tionships between halo structural parameters (e.g., core 
radius and central density). This provides a means to 
evaluate WDM as a solution to the cusp-core problem. 

In SIDM models, the interactions lead to either ther- 
malization or particle loss and the subsequent for- 
mation of a core. Models that have been stud- 
ied are those with large cross sections for sca ttering 
(ISpergel fc Steinhardtil2000t iFirmani et all [2000) or an- 
nihilations (|Kaplinghat et al.ll2000h . In the simplest case 
where the s-wave contributions dominate the cross sec- 
tion, we expect all core densities to cluster around a com- 
mon value. For more complicated models, the core den- 
sity should correlate with the velocity dispersion in the 
core. 

In this Letter we fit high-resolution rotation curves 
of low surface brightness (LSB) disk galaxies using 
cored dark matter density profiles of the type ex- 
pected in WDM and SIDM models. These galaxies 
are dark matter-dominate d down to small radii (e.g., 
Ide Blok fc McGaughl Il996l but see Fuchs 2003), there- 
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fore providing a good laboratory for testing dark matter 
halo predictions. In Section 2, we describe the density 
profiles we use for models of thermal WDM and non- 
thermal WDM from early decays. In Section 3, we fit 
the galaxy data with these halo models and determine 
the sizes of the halo cores. In Section 4, we discuss pre- 
dictions for particle dark matter models based on the 
measured core sizes. A summary is presented in Section 
5. 

2. CORED DARK MATTER HALO MODELS 

Provided below are brief descriptions of the cored halo 
profiles used in our analysis. We consider four differ- 
ent profiles to bracket both a range of behavior in the 
outer region of the density profile and the rapidity of 
the transition from the core to the outer region. In all 
cases, the profiles are set by two independent param- 
eters: a scale den sity and sc ale radius. Recen t WDM 
(|Cohn et al.ll2008l ) and SIDM (|Dave et al.H200l simula- 
tions support this assumption. The specific forms of the 
profiles we use are motivated by arguments in the limit 
of no phase space mixing. However, given the wide range 
of behaviors encapsulated by these profiles, we also use 
them for our analysis of SIDM models. 

2.1. Thermal Warm Dark Matter 

The density profile for the thermal WDM halo, in the 
limit of no phase space mixing, is subject to the con- 
straint that as r — > 0, p — > po[l — (r /r corc ) 2 ]. This en- 
sures that the phase space density is finite for all values 
of total energy This is necessary because the phase space 
density of thermal particles in the early universe is finite 
for all momenta. We assume the following form for the 
thermal WDM halo density profile 



PTh(r) 



po 



(1) 



where po is the central core density and r corc is the core 
radius. 

We test the thermal WDM model with a — 1 and 
a = 2. The a = 1 case can be recognized as the 
cored pseudoisothermal halo traditionally used in rota- 
tion curve fitting. We use this case to test how sensitive 
our results are to the assumed density profile. 

The a = 1 (Thl) and a = 2 (Th2) rotation curves are 



V$ bl (r)=4irGp r, 
V$ h2 (r)=4irGp r, 
with x = rj (y/ar cole ). 
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2.2. Non-thermal Warm Dark Matter from Early 
Decays 

For particles populated by decays of massive particles 
in the early universe, the density profile is subject to the 
constraint that as r — > 0, p — > po(l — r/r corc ) in the limit 
of no phase space mixing. This is a milder core than that 
created by thermal WDM particles. We require that the 



density profile be of the form 



(3) 



where po is the central core density and r corc is the core 
radius. For further discussion about the form of these 
profile s, we refer the reader to Martinez & Kaplinghat 

(pool . 

We test the early decay model with a = 3 and a = 4. 
The choice of these exponents is motivated by the range 
of slopes seen in CDM simulation s in the outer regions 
of the halos ([Navarro et all l2004t iGraham et all 120061 : 
iNavarro et al.ll2008l ). Because the outer regions are built 
up from accreted dark matter particles, there should be 
no differ ence between cold and warm dark matter pre- 
dictions (|Colm et al.ll2008D . 

The rotation curves corresponding to the a — 3 (ED3) 
and a — 4 (ED4) cases are 



^ED3( r )=367rGpor, 

^ED4M=647rGpor c 2 , 
with x = r/(ar corc ). 
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3. CORED DARK MATTER HALO FITS 



We fit the cored dark matter models to the observed 
rotation curves of a sample of nine LSB galaxies. The 
rotation curves are derived f rom high-resolution Ha inte- 
gral fi eld unit spectroscopy (|Kuzio de Narav et al.ll2006l 
2008) and are combined with rotation curves derived 
from l ong-slit Ha spectra (Ide Blok. M cGaugti~ fe Rubinl 
120011: Ide Blok fe Bosnia] I2002f) and Hi ve locity fields 
(|de Blok et al1ll996HStill[l999HSwaters]lll)9l h The inner 
part of the rotation curves set the central core density, 
po, for each halo model, and the data at large radii fix 



PoK 



Because the dark matter is the dominant mass 



component in these galaxies at all radii, we neglect the 
(minimal) contribution to the observed rotation from the 
baryons. To demonstrate that the effect of the baryons 
on the fit parameters is indeed small, we include in the 
figures a representative example (F583-1) where a stel- 
lar mass-to- light ratio based on galaxy color is assumed 
(jKuzio de Narav et al.ll2008[ ). 

In Figure 1 we plot the four cored halo fits over the 
data. The central density and core radius derived for 
each dark matter model are listed in Table 1. 

We find that for most galaxies, the best-fitting halo ro- 
tation curves match the observed galaxy rotation curves 
well. In general, the early decay and thermal WDM mod- 
els produce very similar fits to the data and there is little 
dependence on the value of the exponent a. With the ex- 
ception of the thermal WDM model with a = 2 for some 
of the galaxies, the best-fitting halo rotation curves typ- 
ically overlap at most radii. For the four halo models, 
the sizes of the derived cores are on the order of a kpc 
and span a range of about a factor of 10. 

The observed rotation curves of UGC 4325 and 
DDO 64 are not well-described by the halo models (x^, > 
3). Rather than turning over, the models continue to 
rise past the last observed rotation curve point. The 
implied cores would be larger than the observed radial 
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Fig. 1. — Observed LSB galaxy rotation curves with the best-fitting early decay dark matter (a = 3: solid red, a = 4: dotted orange) 
and thermal WDM (a = 1: short-dash blue; a = 2: long-dash green) halo fits overlaid. (A color version of this figure is available in the 
online journal.) 



TABLE 1 

Best-Fit Cored Halo Parameters and Primordial Phase Space Densities 



Early Decay DM a = 3 Early Decay DM a = 4 



Galaxy 


^core 


PO 


Xr 


Qp 


''"core 


PO 


Xr 




Aftot 


UGC 4325 


4.6 a 


106 b 


3.1 




4.6 a 


106 b 


3.1 






F563-V2 


l.liO.l 


188±30 


0.65 


4.8 


1.3±0.2 


167±25 


0.74 


5.9 


9.8 


F563-1 


1.4±0.1 


106±16 


0.50 


3.2 


1.8±0.1 


90±12 


0.50 


3.1 


14 


DDO 64 


2.7 a 


57 b 


3.3 




2.7 a 


57 b 


3.3 






F568-3 


2.8±0.4 


40±6 


1.5 


0.71 


3.0±0.4 


38±5 


1.4 


1.0 


28 


UGC 5750 


4.3±0.7 


llrfcl 


0.92 


0.42 


4.8±0.8 


10±1 


0.89 


0.53 


30 


NGC 4395 


0.6±0.1 


346±42 


2.6 


21 


0.7±0.1 


278±33 


3.0 


29 


2.6 


F583-4 


0.9±0.1 


98±22 


0.59 


13 


1.1±0.2 


82±18 


0.68 


11 


2.9 


F583-1 


1.7±0.1 


48±4 


0.58 


2.9 


2.1±0.2 


43±3 


0.54 


2.9 


11 



Thermal WDM a = 1 Thermal WDM a = 2 



Galaxy 


7*core 


PO 


Xr 


Qp 


7*core 


PO 


Xr 


Q P 


Aftot 


UGC 4325 


4.1±1.0 


88±6 


3.1 




4.3±1.0 


90±5 


3.1 






F563-V2 


1.5±0.2 


118±18 


0.71 


72 


2.5±0.2 


93±14 


1.4 


14 


4.3 


F563-1 


2.1±0.2 


66±9 


0.43 


35 


3.8±0.2 


47±6 


0.69 


5.9 


7.3 


DDO 64 


2.7 a 


45 b 


3.1 




2.7 a 


47 b 


3.1 






F568-3 


3.8±0.4 


27±3 


1.2 


9.3 


5.0±0.4 


25±2 


1.1 


3.4 


9.2 


UGC 5750 


5.7±0.8 


7.1±0.7 


0.83 


5.4 


7.1±0.7 


7.0±0.2 


0.74 


2.4 


6.9 


NGC 4395 


0.7±0.1 


262±34 


2.9 


478 


1.7±0.1 


121±16 


5.1 


42 


1.7 


F583-4 


1.3±0.2 


66±16 


0.67 


149 


2.4±0.3 


38±8 


1.2 


26 


1.5 


F583-1 


2.5±0.2 


30±2 


0.50 


31 


4.0±0.2 


22±1 


0.77 


7.0 


4.2 



Note. — Best-fit halo parameters (r corc , po), lower limits on the primordial phase 
space densities (Q P ), and the total mass of the system, M to t- The units for r core are 
kpc and the units for po are 10 -3 Mq pc -3 . The units for Q p are 10 -9 M Q pc -3 (km 
s^ 1 )" 3 . The units for Mtot are 1O 1O M . 
a upper limit 
b lower limit 
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Fig. 2. — Mass of each galaxy-halo system as a function of mea- 
sured halo core radius for the Early Decay a = 4 and Thermal 
WDM a = 2 models. Lines of constant Q p for the early decay 



case are overplottcd. 
similar in shape. 



Lines of constant Q p for the thermal case are 
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Fig. 3. — The minimum primordial phase space density, Q p , in- 
ferred from each galaxy in the sample. Each symbol represents a 
different model for the dark matter halo density profile. For a given 
model, these galaxy data are not consistent with a single value of 
Qp. The small gray symbols for F583-1 indicate the results when 
a non-zero stellar mass-to-light ratio is assumed. 

range of the data, and the masses of the systems would be 
about 1O 12 M or larger. Given the poor fits, we exclude 
UGC 4325 and DDO 64 from the remaining analysis. We 
note that the xi values for the NGC 4395 halo fits are 
also large (2.6 - 5.1). These high values are the result 
of the features seen in the rotation curve around 2.5 kpc 
that cannot be fit by smooth dark matter density pro- 
files. However, the velocity profiles do follow the overall 
shape and turnover of the observed rotation curve. 

4. DISCUSSION 

Given a candidate dark matter particle with primor- 
dial phase space density Q p , any halo with total gravi- 
tationally bound mass M to t must have a minimum con- 
stant density core radius that scales inversely with the 



total halo mass as 



-2/3 



M, 



-1/3 



tot 



(Kaplinghat 



120051: iMartinez fc Kaplinghatl 120091 ) . This relation fol- 
lows from dimensional analysis given our assumption 
that the shape of the density profile only depends on 
a scale density and scale radius. It can also be explic- 
itly calc ulated in the c ontext of the "excess mass func - 
tion" of lDehnenl (|2005[ ). as shown bv iKaplinghatl (|2005f ). 
Given the total mass and primordial phase space den- 
sity, a lower bound on the core size can be placed be- 
cause entropy increases or, equivalently, because "phase 
space density" decreases. More correctly, the values of 
the minimum core size reflect a situation where the par- 
ticles have not phase-space mixed. 

The lower limit to the primordial average phase space 
density for thermal WDM models may be written as 



> 



71 



M pc" 3 (km/s)" 



-1/2 -3 
^*core> 



Thl 
Th2 



(5) 



where the upper and lower numbers are for Thl and Th2 
models, respectively, and po is in MqPC -3 and r core is in 
pc. Note that the dependence on the actual shape of the 
halo profile is mild. 

For models of non-thermal WDM from early decays, 
the lower limit to the primordial average phase space 
density may be written in terms of the total mass of the 
dark matter halo as follows 



Q p > 82M Q p C - J (km/s) 



\-3, 



pc 



[ED] (6) 



We calculate the lower limits on Q p for each galaxy and 
each dark matter model using Equations 5 and 6 and list 
them in Table 1. For the ED4 and Th2 models we also 
list the total mass of the galaxy-halo system, M to t, in- 
ferred from the best-fit density profile. The masses deter- 
mined in this way are consistent with and only marginally 
larger than the enclosed mass determined using the last 
observed rotation curve point. The ED3 and Thl models 
formally have divergent masses because of the assumed 
profile shape at large r and are not listed in Table 1. 

The halo parameters presented in Table 1 do not seem 



to obey the r c 



cx M t 



-1/3 



scaling relation that is ex- 



pected if the cores are set by the requirement of no phase- 
space mixing. The observed core radii of our galaxy sam- 
ple span a range of about a factor of ten, which would 
require a factor of ~ 1000 variation in total mass from 
system to system to be explained by the same primor- 
dial phase space density. Given that these galaxies are 
so similar in luminosity and asymptotic rotation speeds 
(Figure 1), this is highly implausible. 

We note that this argument does not rule out the possi- 
bility that the underlying dark matter is warm, but only 
that core sizes are not set by the primordial phase space 
density of dark matter particles. In particular, some of 
the cores could be due to phase-space mixing as a re- 
sult of mergers, but this makes the attribution of cores 
to fundamental dark matter physics ambiguous. In ad- 
dition, we will show that the values of the primordial 
phase space density required are much too small to be 
consistent with constraints arising from the matter power 
spectrum. 

In Figure 2, we present a more detailed comparison 
of theory and data by plotting the inferred total mass 
of each galaxy- halo system, M to t, against the measured 
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Fig. 4. — Halo central density, po, as a function of the maximum 
observed rotation velocity of the galaxy. Each symbol represents 
a different model for the dark matter halo density profile. For a 
given model, po is not constant across the sample, and there is no 
discernible trend in po with V max . The small gray symbols indicate 
the results when a non-zero stellar mass-to-light ratio is assumed. 

halo core radii. For clarity we plot only the results for 
the Early Decay a = 4 and Thermal WDM a = 2 cases. 
Qualitatively similar results are obtained for the ED3 
and Thl cases if we use M v - lr rather than Mtot and do 
not modify the conclusions. For each galaxy we plot the 
combined range of core radii and masses for the two dark 
matter models. For comparison, we have also plotted the 
Mtot vs. (minimum) r core relationship that is expected 
for early-decay dark matter for two choices of Q p . It is 
immediately obvious from Figure 2 that 1) the data span 
a range of only about one order of magnitude in mass, 
2) the data are not consistent with a single value of Q p , 
and 3) mass and core radius are not anti-correlated as 
would be expected from Equation [6l 

The simplest interpretation of this result in the context 
of dark matter models is that the cores in these galaxies 
cannot be set directly by the primordial phase space den- 
sity of dark matter and therefore must be the result of 
baryonic processes. If, however, we insist that a WDM 
model explain these data, then to have a single value 
of Q p for this sample, galaxies with small cores must 
preferentially lose more than 2 orders of magnitude in 
mass, while galaxies with large cores lose very little. This 
is highly unlikely in these undisturbed disk galaxies, as 
feedback from powerful radio sources is observed to occur 
almost a lways only in elliptical galaxies or obvious recent 
mergers (IWilson fc Colbertll995l : riIrrv fc Padovanill995t 
IAntonuccilll993[ ). Additionally, feedback from supernova 
winds is also unlikely to affect these galaxies, as the star 
formation rates in LSBs are known to be lower than 
the rates in high su rface brightness galaxies of similar 
morphological type (iBothun. Impev. fc M cGauah 119971 : 
lO'Neil. Oev. & Bothunl I2007H . We note here that re- 
cent high resolution hydrodynamical simulations have 
produced galaxies with cored CDM halos by includ- 
ing baryonic processes that effectively remo ve mass 
(iGovernato et al.ll2009t iMashchenko et al.ll2008l ). though 
ICeverino fc Klypinl (|2009l ) reach a different conclusion. 



Finally, even if there were a plausible model to explain 
Figure 2, we show below that the required value of Q p is 
in strong disagreement with Lyman alpha forest data. 

In Figure 3, we plot the range of Q p for the galaxies and 
again find that, for a given dark matter model, the data 
are not consistent with a single Q p value. For our sample 
of galaxies, Q p ranges between ~ 10~ 9 and 10 -7 in units 
of M0pc~ 3 (km/s)~ 3 . This result does not change when 
the baryons are accounted for by assuming a non-zero 
stellar mass-to-light ratio, as shown for F583-1 in Figure 
3. These limits on Q p are about 4 orders of magnitude 
smaller than the lower limit on thermal WDM implied 
by the Lyman alp ha forest power sp e ctrum of 10~ 3 
M Q pc- 3 (km/s)- 3 (ISeliak et all [200a IViel~eTaI1 12051 . 
Thus, even if there were a WDM model whose primordial 
phase space density value was in tandem with some other 
process that sets the core sizes in these galaxies, we would 
have a model that is inconsistent with the Lyman alpha 
forest data by orders of magnitude. 

We now consider the SIDM model predictions. This 
is easier to analyze because the SIDM models predict a 
correlation between core size and core density. In most 
models of dark matter with large self-interactions, all 
dark matter halos are predicted to either have the same 
core density or to show a tre nd in pq as a function of ve- 
locity dispersion of the halo dSpergel &: S tcinhardtl l200fl ; 
iFirmani et al.| [2000; Kaplin ghat et al.ll2000l ). One reason 
for this trend is the dependence of the scattering or anni- 
hilation cross section on relative velocity. Additionally, 
adiabatic expansion due to particle loss will result in a 
systematically smaller c ore density in less massive halos 
(jKaplinghat et al.ll2551 . 

We note that a monotonic relation between the cross 
section and the velocity translates to a monotonic rela- 
tion between the core density and the velocity disper- 
sion of the dark matter particles in the core. We expect 
the velocity dispersion in the core to be isotropic and 
proportional to Vmax- It therefore follows that the ex- 
pectation from SIDM models is that the inferred core 
density should be either roughly constant or exhibit a 
monotonic trend with V,nax- We note that if the self- 
interaction process has been operating for differing times 
in these galaxies, for example as the result of a recent ma- 
jor merger, then some dispersion may be introduced into 
the inferred po versus Vmax relation. However, this seem 
unlikely given the uniformity of the sample and the lack 
of observational evidence for any recent disturbance. 

In Figure 4, we plot po against Vm ax for each galaxy 
and show a representative example of how po changes if a 
non-zero stellar mass-to-light ratio is assumed. We find 
that po is not constant across the sample, nor is there 
evidence for a systematic trend in po as a function of 
Vmax- This indicates that the inferred cores in these LSB 
galaxies cannot be directly set by large self-interactions 
(scattering or annihilation) of dark matter. 

5. SUMMARY 

Warm dark matter models and strongly self-interacting 
dark matter models have been proposed to alleviate some 
of the difficulties that CDM faces on small scales. We 
have tested models of thermal WDM and non-thermal 
WDM from early decays with high-resolution rotation 
curves for LSB galaxies. We infer the observed halo core 
radii to span about an order of magnitude around a kpc, 
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while WDM models would predict a spread of only about 
a factor of 2. Additionally, the values of Q p inferred from 
these LSB galaxies are orders of magnitude smaller than 
the lower limits implied by the Lyman alpha forest power 
spectra. Taken together, we interpret these results to 
mean that the cores in these LSB disk galaxies cannot 
be a direct result of WDM particle properties. We also 
find the data to be inconsistent with a single value of core 
density po and find no evidence for a trend in p with the 
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